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ABSTRACT 

Determination of the cosmic infrared background (CIB) at far infrared 
wavelengths using COBEjDIRSE data is limited by the accuracy to which 
foreground interplanetary and Galactic dust emission can be modeled and 
subtracted. Previous determinations of the far infrared CIB (e.g., Hauser 
et al. 1998) were based on the detection of residual isotropic emission in skymaps 
from which the emission from interplanetary dust and the neutral interstellar 
medium were removed. In this paper we use the Wisconsin Ha Mapper (WHAM) 
Northern Sky Survey as a tracer of the ionized medium to examine the effect 
of this foreground component on determination of the CIB. We decompose the 
DIRBE far infrared data for five high Galactic latitude regions into H I™ and Ha- 
correlated components and a residual component. Based on FUSE H 2 absorption 
line observations, the contribution of an ^-correlated component is expected to 
be negligible. We find the Ha-correlated component to be consistent with zero 
for each region, and we find that addition of an Ha- correlated component in 
modeling the foreground emission has negligible effect on derived CIB results. 

Our CIB detections and 2 a upper limits are essentially the same as those derived 
by Hauser et al. and are given by vI v (ESN m -2 sr -1 ) < 75, < 32, 25±8, and 
13±3 at A = 60, 100, 140, and 240 pm, respectively. Our residuals have not been 
subjected to a detailed anisotropy test, so our CIB results do not supersede those 
of Hauser et al. We derive upper limits on the 100 pm emissivity of the ionized 
medium that are typically about 40% of the 100 pm emissivity of the neutral 
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atomic medium. This low value may be caused in part by a lower dust-to-gas 
mass ratio in the ionized medium than in the neutral medium, and in part by a 
shortcoming of using Ha; intensity as a tracer of far infrared emission. If Ha; is not 
a reliable tracer, our analysis would underestimate the emissivity of the ionized 
medium, and both our analysis and the Hauser et al. analysis may slightly over- 
estimate the CIB. We estimate the possible effect for the CIB to be only about 
5%, which is much smaller than the quoted uncertainties. From a comparison of 
the Hauser et al. CIB results with the integrated galaxy brightness from Spitzer 
source counts, we obtain 2 a upper limits on a possible diffuse CIB component 
that are 26 nW m -2 sr -1 at 140 pm and 8.5 nW m~ 2 sr _1 at 240 pm. 

Subject headings: cosmology: observations — diffuse radiation — Galaxy: gen- 
eral — infrared: ISM: continuum — ISM: general 


1. Introduction 

The diffuse cosmic infrared background (CIB) consists of the cumulative energy releases 
in the universe that have either been redshifted, or absorbed and reradiated by dust, into 
the infrared (IR) wavelength region. The CIB therefore provides important constraints 
on the rates of nuclear and gravitational energy release, as well as more exotic forms of 
energy release, over the history of the universe. Over the past few years, analyses of data 
obtained with the Diffuse Infrared Background Experiment (DIRBE) and the Far Infrared 
Absolute Spectrophotometer (FIRAS) onboard the Cosmic Background Explorer ( COBE ) 
satellite have provided the first measurements of the cosmic background in the far infrared 
to submillimeter wavelength region (Puget et al. 1996, Schlegel et al. 1998, Fixsen et al. 
1998, Hauser et al. 1998, Lagache et al. 1999, Lagache et al. 2000). A detailed description of 
the CORE instruments and the CORE mission is given by Boggess et al. (1992), Silverberg 
(1993), Mather, Fixsen, & Shafer (1993), Hauser et al. (1997), and Brodd et al. (1997). 
Recent reviews covering the history of the quest for the CIB, the current detections and 
limits on its spectrum, and the astrophysical implications are given by Hauser and Dwek 
(2001), Kashlinsky (2005), and Lagache, Puget, and Dole (2005). 

The far infrared CIB measurements are limited by the accuracy to which foreground 
interplanetary and Galactic emission can be modeled and subtracted from the COBE data. 
Emission from interplanetary dust is the dominant foreground below about 100 pm and 
emission from interstellar dust is the dominant one at longer wavelengths. 

Different models with different degrees of complexity have been used to remove the 
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interplanetary dust (IPD) emission from the CODE maps. Puget et al. (1996) and Schlegel 
et al. (1998) relied only on the spatial characteristics of the IPD emission and subtracted a 
scaled template based on the DIRBE 25 /xm sky map. The detections reported by Hauser 
et al. (1998) and Fixsen et al. (1998) used the IPD model of Kelsall et al. (1998), which 
was fit to the time variation of the DIRBE data caused by the motion of the earth through 
the interplanetary dust cloud. The uncertainty in the zero level of the emission predicted by 
this model makes a major contribution to the uncertainty of the related CIB measurements. 
Wright (1998) and Gorjian, Wright, and Chary (2000) modeled the IPD using a method 
similar to that of Kelsall et al., with an added constraint that the residual 25 /mi intensity 
after zodiacal light subtraction be zero at high Galactic latitudes. 

Removal of the emission from interstellar dust requires an interstellar medium (ISM) 
template that has a well-defined zero level and correlates well with the spatial variation of 
IR emission in the CODE maps. Before the completion of the Wisconsin H- Alpha Mapper 
(WHAM) survey, observations of Galactic H I emission provided the best ISM template for 
this purpose (Puget et al. 1996; Schlegel et al. 1998; and Hauser et al. 1998). For example, 
Hauser et al. (1998) determined correlations of IR intensity with Galactic H I column density 
for selected regions at high Galactic latitude and high ecliptic latitude, using data in the 
DIRBE 100 /mi, 140 /on, and 240 gm bands. For each region and each wavelength band, the 
H I-correlated component of the infrared emission was subtracted from the data. Careful 
error analysis, including estimates of systematic error in subtraction of the interplanetary 
and Galactic foregrounds, showed that the mean residual intensity was significantly (more 
than 3cr) greater than zero at all three wavelengths. The residual intensity passed tests for 
isotropy at 140 /mi and 240 /mi. so detection of the CIB was claimed at these wavelengths. 

This method of subtracting Galactic foreground emission is subject to error if the ratio 
of Galactic foreground emission to Galactic H I column density varies over the region studied. 
Such variation could occur if there is emission from dust associated with molecular or ionized 
gas, and this emission is not entirely correlated with H I column density. Based on the 
FUSE H 2 absorption line study of Gillmon et al. (2006) and Gillmon and Shull (2006), 
II -2 column density is expected to be negligible compared to H I column density over a 
large fraction of the high latitude sky, and to account for 1% to 30% of total H column 
density for cirrus features brighter than 1.5 — 3 MJy sr -1 in the temperature-corrected 
100 /im map of Schlegel, Finkbeiner, and Davis (1998). Most far infrared CIB determinations 
have used restrictions on Galactic latitude, H I column density, and/or far infrared color 
to exclude lines of sight that may contain significant emission from dust associated with 
molecular gas. On the other hand, the warm ionized medium (WIM) is known to be prevalent 
at high Galactic latitudes. Available data indicates that H II column density is on average 
about one-third of H I column density (Reynolds 1991a), and correlation studies using Ho; as 
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a tracer suggest that a significant fraction of the H II is not spatially correlated with H I (e.g., 
Reynolds et al. 1995, Arendt et al. 1998, Lagache et al. 2000). The ionized gas is expected 
to be subject to approximately the same interstellar radiation field as the neutral gas, and 
the depletion study of Howk and Savage (1999) shows that the dust-to-gas mass ratio in 
the WIM may be similar to that in the warm neutral medium. Thus, the infrared emissivity 
per H nucleus may be similar in the ionized and diffuse H I phases of the ISM, and infrared 
emission from the ionized phase may have a significant effect on any CIB determination that 
is solely based on correlation with an H I template. 

The paucity of tracers of the ionized gas at high latitudes precluded any definitive 
measurement of the amount of IR emission from dust in this gas phase. Nevertheless, sev- 
eral methods were used to estimate and subtract its contribution to the foreground ISM 
emission: (1) Puget et al. (1996) and Lagache et al. (1999) identified the H II emission 
component with a spatially varying component obtained after subtraction of an H I corre- 
lated foreground component from COBE/FIRAS data. This residual component exhibited 
a csc(|6|) dependence (Boulanger et al. 1996), consistent with that expected from a plane 
parallel layer of ionized gas. (2) Fixsen et al. (1998) used maps of H I column density 
and [C II] 158 pm line emission as templates to model the IR emission from the neutral 
and ionized gas phases. They found that essentially none of the high latitude emission ob- 
served by FIRAS correlated with the [C II] template, and their final CIB spectrum was 
more than two times greater than the H II emission subtracted CIB spectrum of Puget et 
al. (1996) at A < 240pm. (3) Arendt et al. (1998) derived an upper limit to the 100 
pm emissivity per H nucleus in the ionized medium from a correlation analysis of Ho:, 
H I, and DIRBE maps for a 10° x 12° region centered at l = 144°, b — —21°. The derived 30- 
upper limit was equal to 3/4 of the 100 pm emissivity per H nucleus for the neutral atomic 
gas in the same region. Assuming this limit is valid for the Lockman hole region, they used 
available Hck and pulsar dispersion measure data for the Lockman hole to place an upper 
limit of 4 nW m -2 sr -1 on the possible contribution of the ionized gas phase to the 100 pm 
foreground emission. This was scaled to obtain upper limits of 5 and 2 nW m~ 2 sr -1 at 
140 and 240 pm, respectively, assuming the spectrum of emission from the ionized phase has 
the same shape as that of the neutral atomic phase. These upper limits axe comparable to 
the overall uncertainties in the CIB determinations. They were noted as possible errors by 
Hauser et al. (1998) but were not included in their quoted CIB uncertainties. If they had 
been included, their reported 240 pm CIB value would still be a 3cr detection. 

The Wisconsin H- Alpha Mapper (WHAM) northern sky survey provided the first Ha 
map sensitive enough to trace the ionized gas phase of the ISM at high Galactic latitude 
(Haffner et al. 2003). Lagache et al. (2000) used prehminary WHAM survey data to 
decompose the 100 to 1000 pm DIRBE and FIRAS data at high latitude into Ha-correlated, 
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H I-correlated, and isotropic components. The regions studied cover about 2% of the sky- 
in the range 25° < |i>| < 50°. They found a significant Ha-correlated component with dust 
temperature of 17.2 K, very similar to that of the H I-correlated component. Assuming a 
constant electron density of 0.08 cm -3 and an electron temperature of 8000 K for the ionized 
gas, they derived an infrared emissivity per H nucleus for this phase that is similar to that 
of the neutral phase. The CIB spectrum they determined from analysis of the FIRAS data 
is consistent with that of Fixsen et al. (1998). The mean residual intensities they found a.t, 
100, 140, and 240 /xm from analysis of the DIRBE data are consistent with the results of 
Hauser et al. (1998), although their uncertainties are larger. 

The agreement between these latest CIB determinations is encouraging, but the dis- 
agreement between the Arendt et al. (1998) and Lagache et al. (2000) results for the 
100 /xm emissivity of the ionized medium is a matter of concern. In this paper, we address 
the possible effects of the ionized medium on the Hauser et al. (1998) CIB results by in- 
cluding an Ha-correlated foreground component in an analysis that is otherwise similar to 
the DIRBE team analysis. We make use of Ha data from the WHAM Northern Sky Survey. 
The paper is organized as follows. In §2 we describe the data sets and sky areas used in 
the analysis. The decomposition method is described in §3. From the decompositions we 
derive the emissivity per H nucleus for the H I and H II phases of the ISM as well as the 
residual emission for the different DIRBE bands. In §4 we present these results and compare 
them with results of previous studies. In §5 we discuss the emissivity results and we discuss 
possible systematic errors if Ha is not a reliable tracer of far infrared emission. Our results 
and conclusions are summarized in 86. 


2. Data Sets 

The regions of the sky analyzed in this paper are shown in Figure 1. Three of them 
were previously analyzed by Hauser et al. (1998) and Arendt et al. (1998): the Lockman 
Hole (LH), a 300 square degree region around the position of lowest H I column density 
at l = 152°, b = +52°; an 8° x 9° region centered on the north ecliptic pole (NEP) at 
l = 96°, b = +30°; and the DIRBE high quality B north (HQBN) region at 6 > +60° 
and f3 > +45°. These regions were originally selected because they were expected to have 
relatively weak Galactic and interplanetary dust foregrounds, and because good quality 
H I observations were available for them. The DIRBE high quality B south region is not 
included in our analysis because it is below the declination limit of the WHAM Northern Sky 
Survey. In addition, we analyze the second quadrant region previously studied by Lagache 
et al. (2000) and a new region in the first quadrant at 30° < l < 80°, 30° < b < 41°. We 
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refer to these as the Q2 region and Ql region, respectively. The Q2 region was originally 
selected because preliminary WHAM Ho; data were available for it. The Ql region was 
selected because it is comparable to the Q2 region in Galactic latitude, it is at ecliptic 
latitude greater than 30 degrees, it has no molecular clouds detected in the CO survey of 
Hartmann, Magnani, and Thaddeus (1998), and it has no cold infrared excess features that 
are characteristic of molecular clouds in the 100 /rm infrared excess map of Reach et al. 
(1998). 

The datasets used in our analysis are listed in Table 1. We use DIRBE 60, 100, 140, 
and 240 /mi mission-averaged skymaps from which the interplanetary foreground emission 
has been subtracted using the model of Kelsall et al. (1998). Foreground Galactic stellar 
emission is negligible at these wavelengths (Hauser et al. 1998, Arendt et al. 1998) and has 
not been subtracted from the data. 

We use H I 21-cm line data integrated over a velocity range that includes all significant 
Galactic emission, converted to H I column density iV(H I) assuming that the line emission 
is optically thin. The H I data for the LH and NEP regions are from Snowden et al. (1994) 
and Elvis, Lockman, and Fassnacht (1994), and were corrected for stray radiation using 
the AT&T Bell Laboratories H I survey (Stark et al. 1992). Estimated 1 <t uncertainties 
in iV(H I) for these regions range from 0.5 x 10 19 cm -2 to 1.0 x 10 19 cm -2 . The H I data 
that we use for the other regions are from the Leiden-Dwingeloo H I survey (Hartmann and 
Burton 1997), which has also been corrected for stray radiation (Hartmann et al. 1996). We 
adopt a lu uncertainty of 1.25 x 10 19 cm -2 for each position in these regions. This is the 
uncertainty of the stray radiation correction that was estimated by Hartmann et al. (1996) 
for a reference position in the Lockman Hole region. 

We use Ha total intensity data from the Wisconsin H- Alpha Mapper (WHAM) North- 
ern Sky Survey (Haffner et al. 2003), which covers the sky north of declination —30°. The 
WHAM instr um ent, has a 1 degree diameter field of view, and the survey was made on a 
regular Galactic coordinate grid with pointings separated by 0.98°/cos b in l and 0.85° in 
b. The spectrum for each pointing was integrated over —80 < vlsr < 80 km s" 1 to obtain 
total Ha intensity. For the regions we study, systematic errors associated with removal of 
geocoronal and atmospheric emission lines from the spectra can be greater than statistical 
measurement uncertainties. These errors can vary from night to night, and sometimes cause 
~ 7° x 7° ’’blocks” of data taken on particular nights to be noticeably offset in mean in- 
tensity relative to their surroundings. We applied offset corrections to affected blocks in 
the HQBN, Lockman Hole, and Q2 regions to remove discontinuities in Ha intensity at the 
block boundaries. For the HQBN and Q2 regions, most blocks appeared to be unaffected 
and these were assumed to set the zero level of the data. For the HQBN region, an offset 
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of 0.15 Rayleigh was added to the data in the areas (102° < l < 117°, 65° < b < 71°), 
(60° < l < 90°, 65° < 6 < 71°), and (l > 117°, 60° < b < 65°). For the Q2 region, 
an offset of 0.3 R was added in (129° < l < 136.6°, 41° < b < 47°) and an offset of 
0.4 R was added in (136° < l < 157°, 47° < b < 51°). For the Lockman Hole, the offset 
corrections are somewhat uncertain and subjective since many blocks appear to be affected, 
by varying amounts. Offset corrections ranged from -0.5 to 0.1 R, and the zero level is deter- 
mined by the WHAM observations of Hausen et al. (2002) for two positions in the region. 
Our decomposition results for the Lockman Hole are consistent with those for the other 
regions, and omitting this region from our analysis does not change our derived GIB results 
significantly. We adopt a la uncertainty of 0.06 Rayleigh for the velocity-integrated, offset- 
corrected Ha intensity at each WHAM pointing. This is near the low end of the range of mis 
Ha dispersion measured within observing blocks at high Galactic latitudes. The dispersion 
tends to increase with increasing mean Ha intensity or decreasing latitude, presumably due 
to increasing rms dispersion in Galactic emission. With the adopted uncertainty, the mean 
Ha signal-to-noise ratio is 4, 8, 18, 18, and 31 for the LH, HQBN, Ql, Q2, and NEP regions, 
respectively. 

The DIRBE data and Leiden H I data were interpolated to the WHAM pointing po- 
sitions for our analysis. The angular resolution of the Elvis et al. and Snowden et al. HI 
data is much better that that of the other datasets, so these data were averaged over the 
WHAM field of view at each WHAM pointing position. Possible contamination by stellar Ha 
absorption was handled either by excluding positions with stars brighter than V = 6.5, or by 
excluding positions with Ha intensity significantly lower than their surroundings. These two 
methods were compared for the LH and NEP regions, and found to give consistent results. 
Some positions with a discrete source detected in the DIRBE data (the planetary nebulae 
NGC 6543, and galaxies NGC 3079, 3310, 3556, 3690, and 4102) were excluded from our 
analysis. The planetary nebula NGC 6210 appears as a bright point source in the WHAM 
data (Reynolds et al. 2005) and was also excluded. 

Maps of the 100 pm intensity, H I column density, and Ha intensity are shown for the Ql 
region in Figure 2, and correlation plots are shown for this region in Figure 3. These figures 
show that the correlation between vl v { 100 pm) and N(H I) is tighter than that between 
n 4,(100 pm) and J(Ha) or that between /(Ha) and N(R I). Similar trends are seen for the 
other regions analyzed in this paper, and for the region around l = 144°, b = —21° studied 
previously (Reynolds et al. 1995, Arendt et al. 1998). 



3. Analysis 


3.1. Decomposition of the Infrared Emission 

The method of analysis follows that used by Arendt et al. (1998) and is similar to that 
used by Lagache et al. (2000). For each DIRBE wavelength band (60, 100, 140, and 240 /mi), 
the infrared intensity distribution within a given region is decomposed into a component that 
is correlated with H I column density, a component that is correlated with Ha intensity, and 
an isotropic component. This is done by making a least squares fit of the form 

vl v { A) = A 1 N{E I) + BJ( Ha) + Cj (1) 

where I v ( A) is the infrared intensity after interplanetary foreground subtraction and Ax, 
E>i, and C\ are fit parameters. Ai is the mean infrared emissivity per H atom for dust 
in the neutral atomic gas phase, B\ is a measure of the mean infrared emissivity of dust 
in the ionized gas phase, and the intercept C\ is the mean residual infrared intensity. For 
comparison, a second decomposition is performed in which an Ha correlated component is 
not included, by making a fit of the form 

uI 1/ (X) = A 2 N(EI)+C 2 . (2) 

Comparison of the derived C\ and values gives the error in the Galactic foreground 
subtraction if the Ha correlated component is neglected. 

Decomposition of an infrared intensity distribution into three components as in equation 
(1) will be successful if the spatial distibutions of iV(H I) and /(Ha) differ significantly from 
each other, and also differ significantly from an isotropic distribution. These conditions are 
met for each of the regions studied here. This is illustrated for the Q1 region by Figures 2 
and 3. The method of analysis also assumes that Ha intensity is a good tracer of far infrared 
emission from the ionized medium. In §5, we discuss possible errors in our results if this is 
not the case. Extinction of the Ha emission is another potential source of error, but its effect 
is negligible for our regions with the data selection criteria described below. We made fits 
to the 100 pm data for each region with and without a correction to /(Ha) for extinction, 
and differences in the results were insignificant. We made the worst-case assumption of pure 
foreground extinction. With the optical depth at Ha calculated as r = 0.04[iV(H I)/10 20 
cm -2 ], the extinction correction was at most 1.22. 

The fits are made using an iterative procedure that minimizes y 2 calculated using mea- 
surement uncertainties in the independent and dependent variables (Press et al. 1992). 
Uncertainties in the fit parameters are determined from the 68% joint confidence region in 
parameter space, using the method of Bard (1974). For each of the regions except for HQBN, 
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data at the highest H I column densities are excluded from the fitting, as described by Arendt 
et al. (1998) for the LH and NEP regions. The zAQlOO //m) - iV(H I) relation deviates from 
linearity in these regions, with excess 100 //in emission relative to N { H I) at the highest H I 
column densities. This type of relation has been found previously for isolated cirrus clouds 
and for large regions of the sky at high Galactic latitude, and the excess 100 //in emission 
has been attributed to emission from dust associated with molecular gas or to non-negligible 
optical depth in the 21-cm line (e.g., Deul and Burton 1993, Reach, Koo, and Heiles 1994, 
Boulanger et al. 1996). For the Lockman Hole, it is consistent with detections of CO line 
emission toward some 100 // m brightness peaks (Heiles, Reach, and Koo 1988; Stacy et al. 
1991; Reach, Koo, and Heiles 1994). We exclude data above H I column densities where 
the ul v ( 1 00 //m) - N(E I) relation begins to deviate from linearity (see Figure 8 of Arendt 
et al. 1998). The cut is made at N(EI) = 1.5 x 10 20 cm -2 , 5.0 x 10 20 cm -2 , 3.0 x 10 20 
cm" 2 , and 3.0 x 10 20 cm -2 for the LH, NEP, Ql, and Q2 regions, respectively. Measurement 
uncertainties are much larger for the DIRBE 140 and 240 fx m bands than for the 60 and 100 
// m bands, so less stringent iV(H I) limits were adopted for the LH and NEP fits in these 
bands, 2.0 x 10 20 cm -2 for the LH and 6.0 x 10 20 cm -2 for the NEP. With these cuts, the 
area of the sky used in the 100 // m analysis is 420, 170, 35, 180, and 380 square degrees for 
the HQBN, LH, NEP, Ql, and Q2 regions, respectively. 


3.2. Limits on Emission from Dust in H 2 

Our analysis does not allow for possible FIR emission from molecular gas that is not 
correlated with H I. This emission is expected to be negligible for the WHAM pointings that 
pass the N(E I) cuts. Gillmon et al. (2006) and.Gillmon and Shull (2006) reported results 
from a FUSE survey of H 2 absorption lines toward 45 AGNs at |6| > 20°. They compared 
their derived values of the molecular fraction, Jh 2 — 2A r (H 2 )/[A ? '(H I) + 2A r (H 2 )], with values 
of temperature-corrected 100 // m intensity D T from the map of Schlegel, Finkbeiner, and 
Davis (1998). (D T is proportional to Galactic dust column density.) The transition from 
low molecular fractions characteristic of optically thin clouds to high values characteristic 
of H 2 self-shielded clouds was found to occur over the range 1.5 < D T < 3 MJy sr -1 , 
with fff 2 varying between 10~ 6 and 10 _1 in this D T range, fn 2 less than 10 -3 at lower D T 
values, and Jh 2 between 10~ 2 and 0.3 at higher D T values. Except for the NEP region, 
most of the WHAM pointings used in our analysis (after the N(H I) cuts have been applied) 
have D T less than 1.5 MJy sr -1 , so fn 2 is expected to be generally less than 10 -3 and the 
uncertainty in derived CIB results due to neglect of this component is negligible. For each 
region except the NEP, we estimate that dust associated with H 2 contributes less than 0.04, 
0.05, and 0.02 nW m -2 sr -1 at 100, 140, and 240 // m, respectively. This assumes that AT( H 2 ) 
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is constant within a region, 2iV(H2) < 10~ 3 x mean iV(H I), and IR emissivity per H nucleus 
in the H 2 phase is given by the slope of the IR - iV(H I) relation. For the NEP, most of 
the WHAM pointings that are used have D T in the transition range from 1.5 to 3.0 MJy 
sr -1 , so fn 2 values as large as 0.1 are possible. However, we find no evidence for significant 
IR emission from H 2 associated dust. The mean residual infrared intensities C\ from our 
analysis of the NEP region are consistent with the values found for the other regions, and 
the 100 pm - N ( H I) relation for the NEP is linear with small scatter over the range of 
iV(H I) used in our analysis (see figures 7 and 8 of Arendt et al. 1998). Excluding the NEP 
region from our analysis would not change our derived CIB results significantly. 


3.3. Test Against Previous Results 

As a check of our analysis method and software, we performed fits of the form of equation 
(1) using data previously analyzed by Lagache et al. (2000) for the Q2 region. Lagache et 
al. performed fits of this form for 122 positions in the region using DIRBE data with 
interplanetary foreground subtracted, Leiden-Dwingeloo H I data, and preliminary WHAM 
data, all smoothed to the ~ 7° resolution of the COBE/ FIRAS. They kindly provided us with 
the data. We performed fits as described above, except to be consistent with the Lagache 
et al. analysis, only measurement uncertainties for the dependent variable (the infrared 
intensity) were used in the calculation of x 2 - Table 2 gives a comparison of our results and 
the Lagache et al. results. The two analyses give the same values for the fit parameters, 
but the values for the fit parameter uncertainties do not agree. Our uncertainties for A\ and 
Bi are larger than those given by Lagache et al. because our uncertainty calculation allows 
for uncertainty due to coupling between the parameters. (The uncertainty quoted for C\ by 
Lagache et al. is not a statistical uncertainty from the fitting, but was determined from the 
distribution of residual intensity values, so comparison with our uncertainty value for C\ is 
not meaningful.) We conclude from this test that differences between our results in §4 and 
those of Lagache et al. are not the result of software errors. 


4. Results 

Parameters from our fits for the five regions are listed in Table 3. Positions with H I 
column density greater that the cuts described in section 3.1 were excluded from the fit- 
ting. For each region and each wavelength, the first row in the table gives results from the 
three-component fit (equation 1) and the second row gives results from the two-component 
fit (equation 2). The uncertainties listed are statistical uncertainties from the 68% joint 
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confidence region in parameter space; the C\ and C 2 uncertainties do not include systematic 
uncertainties that need to be included in determining the total uncertainty for a CIB mea- 
surement. Column 7 of the table lists the number of independent WHAM pointings used for 
each fit. Sample correlation plots showing the fits to the 100 p m data for the Ql and Q2 
regions are shown in Figure 4, and parameters from the three-component fits for each region 
are plotted as a function of wavelength in Figures 5 and 6. 

For all regions in all wavelength bands analyzed, we find that the mean residual intensity 
is nearly the same whether an Ha-correlated component is included in the fitting or not (the 
Ci and C '2 values are in close agreement), and the quality of the fit is nearly the same in the 
two cases. Also, we do not detect significant Ha-con~elated infrared emission; the B x values 
are consistent with zero within the uncertainties, and the A x and A- 2 values are nearly the 
same. The lack of Ha-correlated 100 p m emission is illustrated for the Ql region in Figures 
2 and 3. No correlation is seen between the Ha map in Figure 2c and and the map in Figure 
2d of residual 100 pm emission after subtraction of H I correlated 100 pm emission. This is 
also shown by the correlation plot in Figure 3d. 


4.1. Residual Intensities and CIB Measurements 

Figure 7 shows the residual intensity averaged over the five regions as a function of 
wavelength for each fit type; the C x and C 2 values were each averaged over the regions using 
weighting by 1 /ct 2 . The weighted-average residual intensities for the different fit types agree 
within la statistical uncertainties at all wavelengths, and the agreement is within 2 11 W m~ 2 
sr _1 at 60, 100, and 240 pm. We conclude that addition of an Ha-correlated component in 
modeling the foreground emission at high Galactic latitude has negligible effect on derived 
CIB results. 

To assess whether the weighted-average values for the residual emission can be identified 
as CIB measurements, we calculated the total uncertainty for these values following the 
method used by Hauser et al. (1998). The total uncertainty was calculated as the quadrature 
sum of the statistical uncertainty, interplanetary foreground subtraction uncertainty, and 
DIRBE detector offset uncertainty. Magnitudes of the latter two uncertainties were taken 
from Table 6 of Arendt et al. (1998). (DIRBE gain uncertainty is not included because 
it has the same multiplicative effect on the mean residual and its total uncertainty, and so 
does not affect the signal-to-noise ratio. This uncertainty is 10-14% for the DIRBE bands 
used here.) Based on these total uncertainties, the mean residual emission is more than 3 a 
greater than zero at 140 pm and 240 pm for both the two-component fits and the three- 
component fits, and at 100 pm for the three-component fits. The residual intensity values for 
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the five different regions (Table 3 and Figure 6) are consistent with isotropy at 140 pm and 
240 jjxsx (the values are compatible within their la uncertainties). We have not performed 
detailed anisotropy tests on the maps of residual intensity from our analysis, such as the 
tests performed by Hauser et al. (1998). However, based on the reduced chi-square values 
in Table 3, significant anisotropy is present in the Q1 and Q2 regions at 60 pm and in all 
five regions at 100 /an. 

Table 4 lists the upper limits and detections of the CIB at 60, 100, 140, and 240 /tm from 
Hauser et al. (1998), from our two-component and three-component fits, and from Lagache 
et al. (2000). For cases where the mean residual intensity does not exceed zero by greater 
than 3a or the residual intensity distribution has been shown to be anisotropic, the table lists 
a 2a upper limit on the CIB followed by the mean residual intensity and its la uncertainty 
in parentheses. Our results are consistent with those of Hauser et al. Our residuals have not 
passed a detailed anisotropy test, and in some cases our quoted uncertainties are slightly 
larger than those of Hauser et al., so our CIB results do not supersede the Hauser et al. 
results. 

If Ha intensity is not a good tracer of infrared emission from dust in the WIM, our three- 
component fits may not account for emission from the WIM much better than fits without 
an Ha correlated component do. Any emission from the WIM that is not correlated with 
Ha or H I would contribute to the residual component, so our results would overestimate the 
CIB. The agreement with the Hauser et al. (1998) results could be because their analysis 
overestimates the CIB by a similar amount. In §5.4 we discuss evidence that Ha may not 
be a good tracer, and we estimate the possible effect of emission from the WIM on the CIB 
results from our analysis and from that of Hauser et al. For each analysis, we estimate 
the effect to be only about 5% at 140 and 240 pm, which is much smaller than the quoted 
uncertainties. 


4- 1.1. Systematic Uncertainties 

Although this paper is primarily concerned with the effect of the ionized ISM on CIB 
determination, we discuss here other systematic uncertainties involving the level of the H I 
cut, choice of interplanetary dust model, and choice of photometric calibration. Our results 
for the effect of the ionized ISM have no significant dependence on any of these uncertainties. 

Arendt et al. (1998) found that the sensitivity of the Hauser et al. (1998) CIB results to 
the H I cut is negligible compared to other sources of uncertainty. For the Lockman Hole, for 
example, the mean 100 pm residual intensity was found to vary by less than 0.7 nW m” 2 sr -1 
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when the cut was varied from 1.0 to 2.0 x 10 20 cm -2 . We have investigated the sensitivity of 
the CIB results from our three-component analysis to the H I cut used for our regions that 
were not included in the Hauser et al. analysis, the Q1 and Q2 regions. The effect of varying 
the cut for these regions over the range from 2.5 to 3.5 x 10 20 cm -2 is small. The derived 
mean residual changes by +0.1, -0.7, -0.1, and +0.7 nW m -2 sr _1 at 60, 100, 140, 240 /mi, 
respectively, when the cut is changed from 3.0 to 2.5 x 10 20 cm -2 . These changes are at most 
25% of the uncertainties of the mean residuals given in Table 4. The mean residual changes 
by -0.8, -0.5, -3.0, and -1.3 nW m~ 2 sr -1 in these bands when the cut is changed from 3.0 
to 3.5 x 10 20 cm -2 . These changes are at most 45% of the quoted uncertainties of the mean 
residuals. 

The effect of choice of interplanetary dust model is shown in Table 5. This table com- 
pares mean residual intensities from Hauser et al. (1998) with those obtained by Wright 
(2004). The main difference between these analyses is that Hauser et al. used the interplan- 
etary dust model of Kelsall et al. (1998) and Wright used the model of Gorjian et al. (2000). 
Each of these models was obtained by fitting the time variation observed over the whole 
sky in each of the DIRBE bands with a parameterized model of the dust cloud, but Gorjian 
et al. added a constraint that the residual 25 /an intensity after zodiacal light subtraction 
be zero at high Galactic latitude. This assumed that the 25 CIB is negligible, based 
on CIB upper limits inferred from TeV gamma ray observations of Mrk 501. Table 5 also 
lists the uncertainty of the zodiacal light subtraction for the Hauser et al. analysis, as esti- 
mated by Kelsall et al. This was obtained by comparing results from a series of models with 
different geometries for the density distribution of the dust cloud, which gave comparable 
quality fits to the time variation of the DIRBE data. The difference between the Hauser et 
al. results and the Wright results is not significant relative to the uncertainty in zodiacal 
light subtraction, or relative to the total uncertainty of the mean residual. 

Hauser et al. (1998) noted the effects on CIB results at 140 and 240 /on if the DIRBE 
data are transformed to the FIRAS photometric system. The CIB values are affected by 
differences in zero point and gain between the DIRBE and FIRAS photometric systems, 
which are not significant relative to the zero point uncertainties and gain uncertainties for 
the two systems. Hauser et al. found that transforming the DIRBE CIB results to the 
FIRAS system would reduce the CIB value from 25.0 to 15.0 nW m~ 2 sr -1 at 140 /an, and 
from 13.6 to 12.7 nW rn” 2 sr -1 at 240 /on. The DIRBE team chose to use the DIRBE 
photometric system to report its results since comparing the two photometric systems has 
its own uncertainties associated with the need to integrate the DIRBE map over the FIRAS 
beam, and the FIRAS spectrum over the DIRBE spectral response. The DIRBE team chose 
not to introduce this additional uncertainty, and we made the same choice for this paper. 



-14- 


4-1.2. Comparison with Spitzer Source Counts 

Recently, Dole et al. (2006) used Spitzer MIPS observations to measure the contribution 
of galaxies selected at 24 pm to the CIB at 70 pm and 160 pm. They found that sources 
brighter than 60 p-ly at 24 pm contribute 5.9 ± 0.9 and 10.7 ± 1.6 nW m~ 2 sr -1 to the CIB 
at 70 pm and 160 pm, respectively. With an extrapolation of the 24 jum source counts below 
the 60 pjy detection threshold, they estimated that the full population of 24 pm sources 
contributes 7.1 ± 1.0 and 13.4 ± 1.7 nW m -2 sr -1 at 70 and 160 pm. They noted that these 
should be regarded as lower limits to the CIB since there may be contributions from a diffuse 
background component or from sources missed by the 24 pm selection. They estimated that 
these contributions account for less than ~ 20% of the far infrared CIB. 

A diffuse component of the CIB could be produced by processes such as emission from 
intergalactic dust or radiative decay of primordial particles. Here we compare the Dole et 
al. integrated galaxy brightness at 160 pm with DIRBE CIB measurements at 140 pm and 
240 pm to estimate upper limits on emission from a diffuse component. We scale the quoted 
160 pm intensity of 13.4 ± 1.7 nW m -2 sr -1 to 140 pm and 240 pm using the shape of the 
model spectral energy distribution from figure 13 of Dole et al., which is based on the Lagache 
et al. (2004) galaxy evolution model. This spectrum is also used to apply color corrections 
so the intensities can be compared with the quoted DIRBE CIB values, which assume a 
spectral shape of ul v = constant over the DIRBE bandpass. We obtain integrated galaxy 
intensity values of 13.7 ±1.7 nW m~ 2 sr -1 at 140 pm and 10.7 ± 1.4 nW m~ 2 sr -1 at 240 pm. 
The uncertainties here do not include any uncertainty in the shape of the adopted spectral 
energy distribution. Table 6 compares these values with the DIRBE CIB results, and with 
DIRBE CIB results transformed to the FIRAS photometric system (Hauser et al. 1998). 
The table lists values for the fractional contribution of the integrated galaxy brightness to 
the CIB, and the difference between the CIB and the integrated galaxy brightness. Using 
the CIB results on the DIRBE photometric system yields 2 a upper limits for a diffuse CIB 
component of 26 nW m~ 2 sr -1 at 140 pm and 8.5 nW m~ 2 sr -1 at 240 pm. 


4.2. Infrared Emissivity of the Ionized Medium 

Infrared emissivity results from our three-component fits are shown as a function of 
wavelength in Figure 5. The derived values of emissivity per H atom for the neutral atomic 
gas phase, e(H°), are comparable to previous determinations from high latitude IR - H I 
correlation studies (e.g., Dwek et al. 1997, Reach et al. 1998). The values of emissivity per 
H + ion for the ionized phase, e(H + ), were obtained from the B\ values in Table 3 using a 
conversion factor of /(Ha)/A’(H" f ) = 1.15 Rayleighs/10 20 cm -2 . This is the mean ratio of 
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Ha intensity to pulsar dispersion measure found by Reynolds (1991b) for four high latitude 
lines of sight toward pulsars at z > 4 kpc. It corresponds to an effective electron density, 
n e ff = f n 2 e ds/ f n e d$, of 0.08 cm -3 for an electron temperature of 8000 K and no extinction. 
The J(Ho;)/lV(H + ) ratio ranges from 0.75 to 1.9 Rayleighs/10 20 cm -2 for the lines of sight 
studied by Reynolds, so the uncertainty in the conversion factor is large. The value adopted 
here is the same as that used by Lagache et al. (2000). Figure 5 shows that the derived 
values of e(H + ) are consistent with zero for all regions at all wavelengths. We have checked 
the dependence of derived 100 p m e(H + ) values on the H I cut, and find that they do not 
change significantly when the cut is varied by ± 30%. 

Our derived values of e(H + ) at 100 pm are compared with previous results for other 
regions of the sky in Table 7. All of the results listed are from analyses similar to that used 
here, with either Ha intensity or centimeter wavelength radio continuum intensity used as 
the tracer of the ionized gas. Column 4 lists the electron density adopted in each study for 
calculating the conversion from Ha or radio continuum intensity to N( H + ). Uncertainty 
in this conversion may be as large as a factor of two or more for some regions, but is not 
included in the uncertainty listed for e(H + ). The emissivity values are shown plotted as a 
function of electron density in Figure 8. The 2 a upper limits for our regions are compa.rable 
to that obtained by Arendt et al. (1998) for the 10° x 12° region at l = 144°, b = —21°. 
For the other previously studied regions, 100 pm emission was detected from the ionized 
gas component, and the derived e(H + ) is greater than or equal to the 100 pm e(H°). This 
can be explained if the dust-to-gas mass ratio is the same in the ionized and neutral atomic 
components, with cases of enhanced emissivity in the ionized component caused by Lyman 
alpha heating (for the Barnard’s Loop region, Heiles et al. 2000) or a local source of heating 
(for the Spica region, Boulanger et al. 1995, Zagury, Jones, and Boulanger 1998), or both 
(for the Galactic plane regions, Sodroski et al. 1997). With the exception of the Q2 region, 
the ionized gas in these regions is not representative of the general warm ionized medium 
observed at high latitudes. The electron density is an order of magnitude greater, and the 
ionized gas is much closer to the Galactic midplane (z < 100 pc, compared to an exponential 
z scale height of about 900 pc for the warm ionized medium (Reynolds 1993)). 

The 100 pm e(H + ) value determined by Lagache et al. (2000) for the Q2 region without 
an H I cut is not consistent with our result for the Q2 region with an H I cut applied, or with 
results for the other high latitude regions that sample the general warm ionized medium. 
We have applied our three-component decomposition to the full Q2 region (with no H I cut) 
at 100 pm at the 1 degree resolution of the WHAM data. This gave e(H + ) = 8.5 ± 1.0 
nW m -2 sr _1 /10 20 cm -2 , which is comparable to the upper limits in figure 8, and e(H°) = 
18.7 ± 0.4 nW m -2 sr _1 /10 2 ° cm -2 . (If an extinction correction is applied to the Ha data 
as described in §3.1, the derived e(H + ) is 5.5 ± 1.0 nW m -2 sr -1 /10 20 cm -2 .) The difference 
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between these results and those of Lagache et al. is probably due to the different angular 
resolution used. Their analysis used data for 122 positions at 7 degree FIRAS resolution, 
but the data were oversampled and there are only about 30 independent FIRAS pointings 
within the region. Our analysis used data for 1762 independent WHAM pointings that are 
within 3.5 degrees (half of the FIRAS beam width) of any of their 122 positions. It is possible 
that some unmodeled effect, such as dust associated with molecular gas, optical depth in the 
21-cm line, or variation of e(H°), happens to correlate with the Ha data at FIRAS resolution, 
but does not correlate as well at WHAM resolution. The u I u (100 pan) - A(H I) relation for 
the region shows curvature, which suggests that molecular gas or 21-cm line opacity may 
be present. We find that 40% of the region has temperature-corrected 100 pan intensity 
greater than the 3 MJy sr -1 threshold for significant fractional H 2 abundance (Gillmon and 
Shull 2006). For these reasons, we consider the emissivity per H + ion derived for the full Q2 
region, from either our analysis or Lagache et al.’s analysis, to be of questionable reliability. 


5. Discussion 

Our derived emissivity values for the ionized medium are statistically consistent with 
zero. From Table 7, our derived 2a upper limits on 100 pan emissivity per H + ion are 
0.33, 1.11, 0.37, 0.57, and 0.30 of the 100 pan emissivity per H atom for the HQBN, LH, 
NEP, Ql, and Q2 regions, respectively. We adopt 0.4 as a representative upper limit on 
e(H + )/e(H°) at 100 pan for these regions. Possible explanations for this low value include a 
lower dust-to-gas inass ratio or a weaker radiation field in the ionized medium than in the 
neutral medium, a difference in grain optical properties or grain size distribution, an error 
in our adopted /(Ha)/iV(H + ) conversion factor, or a shortcoming of using Ha as a tracer of 
infrared emission from dust in the ionized medium. In this section, we discuss some of these 
possibilities. Based on available observations and models, it appears that our low derived 
e(H + )/e(H°) ratio is partly due to a lower dust-to-gas mass ratio in the WIM and partly due 
to error in using Ha as a tracer, but one or more of the other factors may also contribute. 
We discuss possible implications for derived CIB results in §5.4. 


5.1. Dust-to-gas Mass Ratio 

A lower dust-to-gas mass ratio would be expected if the ionized gas in these regions has 
greater z extent and lower density than the neutral gas, as is typical for the WIM in the 
solar vicinity (Reynolds 1991a). From interstellar absorption line observations, it has been 
inferred that abundances of heavy elements in the form of dust decrease with increasing 
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2 , and also decrease with decreasing gas density (e.g., Savage and Sembach 1996). Most 
of these studies have pertained to the neutral atomic medium, but from observations of 
A1 III and S III absorption lines, Howk and Savage (1999). found evidence that about 60- 
70% of aluminum atoms in the WIM are in dust, compared to about 90% of A1 in dust in 
H II regions at low 2 . This result is based on observations of two lines of sight that sample 
the WIM up to z distances of 690 pc and 2800 pc, and four lines of sight through low density 
(n e ~ 0.2 to 4 cm -3 ) H II regions at z < 200 pc. Howk and Savage also noted that the 
dust phase A1 abundance they obtained for the low z H II regions is comparable to previous, 
somewhat uncertain determinations for the warm neutral medium at low 2 . Assuming that 
this low 2 dust phase A1 abundance is valid for the neutral atomic medium in the regions 
studied here, and that 100 /xm emissivity per H nucleus varies in proportion to A1 dust 
phase abundance, one would expect e(H + ) to be about 20-30% lower than e(H°) at 100 /xm. 
Thus the Howk and Savage results suggest that the difference in dust abundance between the 
neutral and ionized components is not great enough to fully explain our derived upper limit 
on e(H + )/e(H°). To confirm this,. dust phase A1 abundance determinations for additional 
lines of sight through the WIM would be of interest, as would calculations of the relation 
between A1 dust phase abundance and 100 /xm emissivity for dust grain models. 


5.2. Interstellar Radiation Field 

The lower 100 /xm emissivity derived for the ionized medium could also be explained 
if it were subject to a. weaker interstellar radiation field. Heating of dust by Lyman alpha 
is expected to be negligible at the electron density estimated for the WIM (Spitzer 1978, 
Heiles et al. 2000). Also, there are no early type stars that are close enough to the lines 
of sight through our regions to cause significant dust heating relative to that of the general 
interstellar radiation field. Models of the spatial distribution of the interstellar radiation 
field at visual and ultraviolet wavelengths predict that the mean intensity increases with 
increasing z up to about 200 pc, as light from distant stars in the Galactic disk becomes 
less attenuated, and then decreases with further increase in 2 due to geometrical dilution. 
Wakker and Boulanger (1986) used their model of the radiation field to calculate the expected 
100 /xm intensity of a diffuse cloud with a mixture of silicate and graphite grains. The 
intensity was calculated for different distances of the cloud along two high latitude lines 
of sight (toward b = 90°, and towai'd l =■ 180°, b = 60°). They found that the 100 /xm 
intensity varies by less than ±20% over the range in cloud height from 2 = 0-1 kpc. Thus, 
it appears that the radiation field does not vary enough to explain our derived limit on 
e(H + )/e(H°), even if all the H I gas were located near the maximum of the radiation field 
and all of the H II gas were at 2 = 1 kpc. 
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5.3. /(Ha) /IV(H + ) Conversion Factor 

To change our derived upper limit on e(H + )/e(H°) at 100 //m from 0.4 to 1.0, the adopted 

1 (Ha) / N (R + ) conversion factor would need to be changed from 1.15 to 2.9 Rayleighs/ 10 20 

cm -2 . This is significantly larger than the largest value of 

1.9 Rayleighs/ 10 20 cm -2 measured by Reynolds (1991b) for lines of sight to four pulsars 
at high 2 . Values obtained by Arendt et al. (1998) for an additional 5 pulsar lines of sight 
range from 0.8 to 1.2 Rayleighs/ 10 20 cm -2 . Thus, it appears that there is not an error in 
the adopted conversion factor value that is large enough to fully explain the low emissivity 
derived for the ionized medium. 

Independent evidence that supports our adopted conversion factor comes from interstel- 
lar absorption line studies. Our adopted value corresponds to an effective electron density 
of 0.08 cm -3 . Electron density estimates from observations of absorption lines of excited C + 
toward extragalactic objects and high 2 stars are comparable to this. In the most extensive 
study to date, Lehner, Wakker, and Savage (2004) presented results for 43 such lines of sight 
at | 6 1 > 30°. Most of the observed absorption line components are at low velocity. For these, 
they find a mean density of (n e ) = 0.08 ± 0.04 cm -3 (la dispersion). For the Intermediate 
Velocity Arch, they find (n e ) = 0.03 ±0.01 cm -3 , probably lower because the gas is at higher 

2 (<■'-' 1 kpc). The derived n e values are averages over C + regions in both the warm ionized 
medium and the warm neutral medium, but Lehner et al. used WHAM Ha data to estimate 
that at least 50% of the excited C + column density originates in the WIM for an average 
line of sight. 


5.4. Ha as a Tracer of Infrared Emission 

Our derived emissivity values and CIB results are subject to error if the WHAM Ha 
data are not a good tracer of far infrared emission from the WIM in the regions we study. 
Far infrared intensity is proportional to dust column density, whereas Ha intensity is pro- 
portional to the square of the ionized gas density integrated along the line of sight, /(Ha) 
oc T e ~°- 92 f n 2 c/s, where n e is electron density and T e is electron temperature (Reynolds 1992). 
Thus errors would be expected in our results if the spatial variation of Ha intensity for a 
region is caused more by differences in mean electron density or mean electron temperature 
for different lines of sight than by differences in ionized gas column density. The approximate 
esc \b\ dependence of high latitude WHAM data (Haffner et al. 2003) provides evidence that 
Ha is a reasonable tracer of A r (H + ) on large angular scales, but this isn’t necessarily true on 
the scales within the regions studied here. 
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Haffner, Reynolds, and Tufte (1999) and Reynolds, Haffner, and Tufte (1999) have found 
evidence that variations in Ha intensity may be largely due to variations in electron density, 
based on their observations of Ha, [N II] A 6583, and [S II] A 6716 line intensities in the 
region 123° < l < 164°, —35° < b < —6°, and previous observations of these lines in halos of 
edge-on galaxies. The [N II] /Ha and [S II] /Ha intensity ratios are observed to increase with 
increasing z, while the [S II] /[N II] ratio is nearly constant. Interpreting the variations in [N 
II] /Ha as primarily due to variations in electron temperature, Haffner et al. and Reynolds 
et al. inferred T e values ranging from 6000 K to 11000 K, with temperature increasing as Ha 
intensity decreases. They showed that this anticorrelation can be explained if variations in 
Ha are largely due to variations in mean electron density, and a supplemental source of gas 
heating is present that dominates over photoionization at low density, causing temperature 
to increase with decreasing density. A number of possible supplemental heating mechanisms 
have been proposed, and Reynolds et al. estimate the heating rate that would be needed to 
explain the observations for each mechanism. 


If Ha is not a good tracer of ionized gas column density, our method of analysis tends 
to underestimate the infrared emissivity per H + ion and overestimate the CIB. To show this, 
we consider a simple model in which the Ha variation in a region is partly due to variation 
of N( H + ) and partly due to variation of effective electron density n e //, 

/(Ha) jV(H+) n ,„ 

(/(Ha)) (iV(H+)> (n,„) W 


where. (I (Ha)), (N( H + )), and (n e //) are averages over all lines of sight through the region. 
This assumes that electron temperature is constant, and (N( H + ) n e ff) — (A r (H + )) ( n e ff ). 
We adopt 


N( H+) 
(JV(H+)) 


I (Ha) 
(/(Ha)) J 


( 4 ) 


and 


n eff 


/(Ha) 


neff) L(/(Ha))J 


1 "V 


( 5 ) 


where p is the fraction of the variation of log / (Ha) that is caused by variation of N (H + ) . We 
assume that, averaged over large angular scales, N( H + ) and /(Ha) are directly proportional, 


(N( H + )) = <I(Ha))/c, (6) 

where c is our adopted conversion factor of 1.15 Rayleighs/ 10 20 cm -2 . The infrared emissivity 
per H nucleus is assumed to be constant within each gas phase, so the infrared emission from 
each phase is proportional to its gas column density. The infrared emission from the ionized 
phase is 

v iwm = e ( H +)jV(H + ) = e(H + )(/(Ha)) 1 - p /(Ha)Vc. (7) 
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Assuming the distributions of N( H I) and /(Ha) are uncorrelated, the Ho; coefficient B x 
that would be obtained from our decomposition is the mean slope of the uI^ IM - /(Ha) 
relation. If the distribution of Ha intensities is symmetric about the mean, this is given by 


B\ = 


dvl ™ 
dl{ Ha) 


= pe(H + )/c 

/(Ha) = (I (Ha)) 


( 8 ) 


and the emissivity per H + ion derived from the analysis underestimates the actual emissivity, 


^<teru)ed(hf ) cBx pe(H ). 


( 9 ) 


Thus, our derived limit e(H + )/e(H°) <0.4 may be consistent with no real difference in 
emissivity between the ionized and neutral phases if only a small fraction of the Ha variation 
in our regions is caused by variation of iV(H + ), i.e., if p < 0.4 in the context of this simple 
model. Berkhuijsen, Mitra, and Muller (2006) have estimated p = 0.68 ± 0.04 for the high 
latitude diffuse ionized gas, using pulsar dispersion measure data and WHAM Ha data 
toward a sample of 157 pulsars at \b\ > 5° (see their Figure 7b). This result suggests that 
our low derived e(H + )/e(H°) is at least partly due to error in using Ha as a tracer. 

The overestimate of the CIB for our simple model is given by the intercept of a linear 
fit to the vlf 1M - /(Ha) relation. To a good approximation, this intercept is given by 

A(^) = - B,(/( Ha)). (10) 

Using equations (6), (7), and (8), we obtain 


A (1/4,) = (1-p)£(H+)(1V(H + )). (11) 

We have used this result to estimate the possible effect of emission from the WIM that is not 
correlated with Ha or 1V(H I) on the CIB results derived from our three-component fits. We 
subtracted A (vl u ) offsets from the derived residual intensities C x for each region assuming 
p = 0.5±0.5, (jV( H+)> = (/(Ha))/(1.15 Rayleighs/10 20 cm" 2 ), and e(H+) = (0.5A0.5) e(H°), 
using values of A x in Table 3 for e(H°). For each wavelength, we then calculated the weighted 
average of the reduced residual intensity over the five regions, treating the ±0.5 uncertainties 
in p and e(H + )/e(H°) as additional independent sources of error. The resulting average 
residual intensities are 21 ±27, 17.0 ±6.3, 23.5 ±8.5, and 12.0 ±2.9 nW m~ 2 sr _1 at 60, 100, 
140, and 240 pm, respectively. At 140 and 240 pm, these results are only 6% lower than the 
results from our three-component fits. 

Equation (11) is the same as the expression for the CIB overestimate from a two- 
component fit (using H I without Ha) where p is the fraction of iV(H + ) that is correlated 
with N( H I). Thus we can use the same procedure to estimate the possible effect of emission 
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from the WIM that is not correlated with iV(H I) on the CIB results of Hauser et al. (1998). 
We make the same assumptions as in the previous paragraph for p, (N( H + )) f and e(H + ). 
The HQBS region included in the Hauser et al. analysis is below the declination limit of the 
WHAM survey, so for this region we used data from the Southern H-Alpha Sky Survey of 
Gaustad et al. (2001) as processed by Finkbeiner (2003). We find that the weighted-average 
residual intensities from the Hauser et al. analysis would be reduced to 20.9 ±6.2, 23.6 ±7.0, 
and 12.9 ± 2.5 nW m -2 sr -1 at 100, 140, and 240 pm, respectively. These results are only 
about 5% lower than those of Hauser et al. 

The possible effect of the ionized medium is also estimated to be small for the 125- 
2000 pm CIB spectrum determined by Fixsen et al. (1998) from FIRAS observations. In 
one of their analyses, Galactic emission template maps were constructed from DIRBE 140 pm 
and 240 pm data with the Hauser et al. (1998) CIB and zodiacal light subtracted, and the 
CIB spectrum was obtained by correlating FIRAS data with these templates. Our estimate 
of ~ 5% for the possible error of the DIRBE CIB results also applies to the 140 A t 240 pm 
part of the CIB spectrum derived from this analysis. At longer wavelengths, the error is 
expected to decrease. This is because the spectrum of emission from the ionized medium 
is expected to be similar to the spectrum of emission from the H I phase, and the ratio of 
this spectrum to the CIB spectrum decreases with increasing wavelength for A > 240 pm. 
Since all three of the Fixsen et al. analyses gave consistent results for the CIB spectrum, 
the possible error due to the WIM should be small for their average CIB spectrum from the 
three methods. 


6. Summary 

We used WHAM Hu data as a tracer of far infrared emission from the warm ionized 
phase of the ISM in an effort to determine the intensity of this emission at high Galactic 
latitudes and to assess its effect on determination of the cosmic far infrared background. 
We studied five high latitude regions, including regions previously analyzed by Hauser et 
al. (1998) and Lagache et al. (2000). For each region, we decomposed COBE/DIRBE data 
at 60, 100, 140, and 240 pm into a sum of an H I correlated component, an Ha correlated 
component, and a residual component. Uncertainties in our results due to omission of an H 2 
correlated component are expected to be negligible, based on results of a F USE high latitude 
H 2 absorption line survey. We found that the intensity of the Ha correlated component is 
consistent with zero within the uncertainties for all regions at all wavelengths. From the 
mean intensities of the residual components, we derived estimates of the CIB at 140 and 
240 pm and upper limits to the CIB at 60 and 100 pm (Table 4). We repeated the analysis 
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without including an Ha correlated component, and the derived CIB results did not change 
significantly. Our CIB estimates and upper limits are similar to previous CIB determinations 
for which the FIR emission from the ISM was traced only by H I column density. We conclude 
that addition of an Ha correlated component in modeling the ISM emission at high Galactic 
latitude has negligible effect on derived CIB results. We did not perform detailed anisotropy 
tests on the maps of residual intensity from our analysis, so our CIB results do not supersede 
the results of Hauser et al. (1998). 

We derived 2a upper limits to the 100 pm emissivity per H + ion for the five regions 
that are typically about 40% of the emissivity per H atom for the neutral atomic medium. 
Available evidence suggests that this low value is partly due to a lower dust-to-gas mass ratio 
in the ionized medium than in the neutral atomic medium, and partly due to a shortcoming 
of using Ha as a tracer of FIR emission, which causes our analysis to underestimate the 
emissivity of the ionized medium. Other possible effects that may play a role include a 
weaker radiation field in the ionized medium than in the neutral medium, a difference in 
grain optical properties or grain size distribution, or an error in our adopted /(Ha)/iV(H + ) 
conversion factor. (The value of 100 pm emissivity per H + ion derived by Lagache et al. 
(2000) is much greater than the upper limit we derived for their region. Our analysis differs 
from theirs in that (1) we exclude positions where iV(H I) is greater than 3 x 10 20 cm -2 , 
where emission from H 2 associated dust or 21-cm line opacity may be significant., and (2) we 
analyze data at a resolution of 1° instead of T°.) 

Ha observations have previously been used in this kind of analysis with apparent success 
for ionized regions that have higher density and are at low £. However, for the general high 
latitude WIM, evidence from Reynolds, Haffner, and Tufte (1999) and Berkhuijsen, Mitra, 
and Muller (2006) suggests that variations in Ha are not entirely due to variations in H + 
column density, but are also due to differences in mean electron density for different fines of 
sight. Thus, Ha may not be an accurate tracer of far infrared emission from the WIM. 

If Ha intensity is not a good tracer, any emission from the WIM that is not correlated 
with Ha or H I would contribute to our derived residual component for each region, so 
our analysis could overestimate the CIB. The agreement with the Hauser et al. (1998) 
results could be because their analysis overestimates the CIB by a similar amount. We used 
WHAM data to estimate the possible effect on our CIB results and on the Hauser et al. 
results, assuming that the mean Ha intensity for each region can be used to estimate its 
mean H + column density. In each case, we estimated the effect to be only about 5% at 
140 and 240 pm, which is much smaller than the quoted uncertainties. The possible effect 
of emission from the WIM is also estimated to be small for the 125-2000 pm CIB spectrum 
determined by Fixsen et al. (1998) from FIRAS observations. 
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We estimated upper limits on a possible diffuse component of the CIB by comparing 
the Hauser et al. (1998) CIB results with the integrated galaxy brightness determined by 
Dole et al. (2006) from Spitzer source counts. We obtained 2 a upper limits on a diffuse 
component of 26 nW m~ 2 sr _1 at 140 pm. and 8.5 nW m~ 2 sr _1 at 240 pm. 

We thank F. Boulanger, J. C. Howk, and G. Lagache for useful discussions, and G. 
Lagache for providing us with data from the previous study of the Q2 region. We thank 
the referee for helpful comments. This paper is dedicated to the memory of our friend and 
colleague Thomas J. Sodroski, and we acknowledge his contributions in the formative stages 
of this work. This research was supported by NASA Astrophysical Data Program NRA 
99-01-ADP-137. 
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Table 1. Datasets 


Dataset 

Resolution 

Velocity 

Integration 

Reference 

DIRBE Zodi-Subtracted 
Mission Average Maps 

0?7 


Hauser et al. (1997) 

Leiden— Dwingeloo 
H I Survey 

0?6 

—450 < ulsr < 400 km s" 1 

Hartmann Sz Burton (1997) 

Lockm&n Hole 
H I map 

21' 

—150 < vlsr < 100 km s _1 

Snowden et al. (1994) 

NEP H I Map 

21' 

—150 < ulsr < 150 km s" 1 

Elvis et al. (1994) 

WHAM Ha 
Sky Survey 

1?0 

—80 < vlsr < SO km s“ x 

Haffner et al. (2003) 


Table 2. Test of Analysis for Q2 Region 


Wavelength 

(H 

Ai a 

(nW m -2 sr~ 1 /10 20 cm -2 ) 

Hi 

(nW m~ 2 sr~ x / Rayleigh) b 

Ci 

(nW m~ 2 sr _1 ) 

Source 

100 

14.9 ±0.1 

14.0 ±0.4 

23.4 ±6.3 

Lagache et al. (2000) 

100 

14.9 ±1.1 

14.0 ± 3.0 

23.4 ±3.8 

This paper 

140 

20.0 ±0.6 

23.4 ± 1.9 

24.3 ± 11.6 

Lagache et al. (2000) 

140 

20.0 ±2.0 

23.4 ± 5.9 

24.3 ± 7.9 

This paper 

240 

9.6 ±0.3 

12.7 ±0.8 

11.0 ±6.9 

Lagache et al. (2000) 

240 

9.6 ±1.3 

12.7 ±3.8 

11.0 ±5.1 

This paper 


Parameter values from fits of the form i/I v (X) = AiA r (HI) + Bi/(Ha) + C\ 
b l Rayleigh = 10 s /47 r photons s" 1 cm“ 2 sr” 1 = 0.24 nW m“ 2 sr^ 1 at Ha 
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Table 3. Decomposition Results 


Region 

DIRBE Band 

Qxm) 

Fit type 

i 

(nW ra“ 2 sr” 1 /10 20 cm” 2 ) 

Bi 

(nW m"“ 2 sr” 1 / R.ayleigh) c 

(nW rn“ 2 sr -1 ) 

■^points 

xl 

HQ BN 

60 

l a 

8.4 ± 0.6 

0.5 ± 2.3 

18.0 ± 1.0 

535 

1.25 



2 b 

8.4 db 0.5 


18.2 ± 0.6 

535 

1.25 


100 

1 

23.1 ± 1.2 

-5.3 ± 5.9 

16.8 ± 2.3 

535 

2.11 



2 

22.7 ± 0.9 


14.8 ± 1.3 

535 

2.12 


140 

1 

19.6 ± 7.1 

-30.8 ±30.0 

36.0 ± 12.9 

535 

1.31 



2 

17.3 ± 5.5 


23.8 ± 7.7 

535 

1.32 


240 

1 

8.4 ± 2.1 

1.8 ± 8.5 

11.2 ± 3.6 

535 

1.12 



2 

8.6 ± 1.7 


11.9 ± 2.3 

535 

1.12 

LH 

60 

1 

4.9 ± 1.1 

4.8 ± 2.4 

24.2 ± 1.0 

215 

0.86 



2 

5.3 ± 0.9 


25.0 ± 0.9 

215. 

0.91 


100 

1 

18.8 ±2.0 

7.7 ± 5.1 

19.3 ±2.0 

215 

2.20 



2 

20.2 ± 1.7 


19.9 ± 1.7 

215 

2.27 


140 

1 

18.6 ± 11.8 

-8.6 ± 39.6 

24,6 ± 13.1 

246 

1.34 



2 

17.7 ± 9.4 


23.2 ± 10.6 

• 246 

1.34 


240 

1 

8.9 ± 3.4 

-18.0 ± 10.6 

17.7 ± 3.8 

246 

1.02 



2 

7.3 ± 2.8 


14.7 ± 3.2 

246 

1.05 

NEP 

60 

1 

8.7 ± 0.8 

-0.1 ± 1.2 

14.1 ± 2.8 

45 

1.12 



2 

8.7 ± 0.6 


14.2 ± 2.2 

45 

1.10 


100 

1 

18.9 ±2.0 

0.7 ± 2.7 

16.4 ± 6.6 

45 

2.86 



2 

19.3 ± 1.3 


16.0 ± 5.2 

45 

2.81 


140 

1 

8.9 ± S.2 

9.0 ± 11.5 

54.8 ± 29.8 

54 

1.02 



2 

13.6 ± 5.5 


50.5 ±23.9 

54 

1.04 


240 

1 

4.6 ± 2.8 

5.8 ± 4.2 

19.9 ± 10.3 

54 

1.10 



2 

7.5 ± 2.0 


17.5 ± 8.6 

54 

1.22 

Qi 

60 

1 

6.6 ± 0.8 

1.7 ± 1.4 

17.8 ± 2.0 

217 

2.51 



2 

6.9 ± 0.6 


19.1 ± 1.3 

217 

2.53 


100 

1 

20.8 ± 1.7 

2.8 ± 3.7 

13.9 ± 4.6 

225 

4.22 



2 

21.1 ± 1.3 


16.0 ± 2.8 

225 

4.23 


140 

1 

24.1 ± 7.0 

-13.3 ± 14.0 

38.0 ± 18.6 

234 

1.33 



2 

22.4 ± 5.5 


27.7 ± 11.9 

234 

1.34 


240 

1 

11.0 ± 2.3 

2.6 ± 4.5 

8.5 ±6.3 

235 

1.42 



2 

11.3 ± 1.9 


10.5 ± 4.0 

235 

1.42 

Q2 

60 

1 

5.4 ±0.5 

0.3 ± 0.8 

22.5 ± 1.1 

499 

2.61 



2 

5.5 ± 0.4 


22.6 ± 0.8 

499 

2.60 


100 

1 

15.4 ± 0.9 

-1.0 ± 2.5 

21.9 ± 1.4 

487 

6.63 



2 

15.2 ± 0.6 


21.6 ± 1'.4 

487 

6.63 


140 

1 

24.5 ± 4.8 

-0.9 ± 9.4 

15.3 ± 8.6 

496 

1.21 



2 

24.3 ± 3.4 


14.9 ± 7.4 

496 

1.21 


240 

1 

10.2 ± 1.6 

-0.4 ± 3.1 

11.1 ± 3.1 

497 

1.33 



2 

10.1 ± 1.2 


10.9 ± 2.6 

497 

1.33 


a Ftt type 1: vl v { X) - AiiV(HI) + BiJ(Ha) + C\ 
b Pit type 2: vl v { X) = ^ 2 ^(HT) + C 2 

C 1 Rayleigh = 10 6 /4tp photons s"" 1 cm -2 sr -1 = 0.24 nW m“ 2 sr — 1 at Ha 
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Table 4. CIB Limits and Detections 21 


ISM 


i/7^ (CIB) (nW m -2 sr" 1 ) 



Tracers 

60 fxm 

100 fim 

140 fim 

240 fi m 

Reference 

H I 
H I 

H I & Ha 
H I & Ha 

< 75 (20.6 ± 27) 

< 75 (20.5 ± 27) 

< 75 (20.9 ± 27) 

< 34 (21.9 ±6.1) 

< 30 (18.2 ±6.1) 

< 32 (19.8 ±6.1) 
23.4 ± 6.3 b 

25.0 ±6.9 

22.0 ± 7.0 6 

25.1 ±8.0 6 

< 47 (24.2 ± 11.6) 

13.6 ±2.5 
12.1 ±2.5 6 
12.8 ± 2.8 b 
< 25(11.0 ±6.9) 

Hauser et ah (1998) 
This paper 
This paper 
Lagache et ah (2000) 


a Upper limits are given for cases where the residuals are not 3 a greater than zero or where they are non-isotropic. 
b A test for anisotropy was not performed. 


Table 5. Effect of Interplanetary Dust Model on CIB Determination 




i //„ (nW 

m 2 sr -1 ) 


60 fim 

100 ptm 

140 /7m 

240 /7,m 

Hauser et ah (1998) mean residual a 

20.6 ± 27 

21.9 ± 6.1 

25.0 ± 6.9 

13.6 ± 2.5 

Wright (2004) mean residual b 

-8 ± 14 

12.5 ± 5 

22 ± 7 

13 ± 2.5 

Residual difference (row 1 - row 2) 

28.6 

9.4 

3 

0.6 

Zodi subtraction uncertainty c 

26.7 

6.0 

2.3 

0.5 


a Based on the IPD model of Kelsall et al. (1998). 
b Based on the IPD model of Gorjian et ah (2000). 

C lcr uncertainty in zodiacal light subtraction for the Hauser et ah analysis, from Kelsall et 
ah (1998). 


Table 6. Comparison of Integrated Galaxy Brightness and CIB Measurements 


140 240 fj , m 


Integrated galaxy brightness a b 
DIRBE CIB a c 

DIRBE CIB transformed to FIRAS scale a c 
Integrated galaxy brightness/CIB (DIRBE scale) 
Integrated galaxy brightness/CIB (FIRAS scale) 
CIB (DIRBE scale) - Integrated galaxy brightness a 
CIB (FIRAS scale) - Integrated galaxy brightness a 


13.7 ± 1.7 

25.0 ± 6.9 

15.0 ± 5.9 
0.55 ± 0.17 
0.91 ± 0.37 

<26 (11.3 ±7.1) 
< 14 (1.3 ±6.1) 


10.7 ± 1.4 

13.6 ± 2.5 

12.7 ± 1.6 
0.79 ± 0.18 
0.84 ± 0.15 

< 8.5 (2.9 ±2.8) 

< 6.2 ( 2.0 ± 2 . 1 ) 


a yJt/ in nW m““ 2 sr” 1 . 

b Dole et ah (2006) 160 jam result from Spitzer MIPS observations, scaled to DIRBE wave- 
lengths. 


; From Hauser et ah (1998). 
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Table 7. 100 fxm Emissivity Determinations 


Region 

1 

M 

Adopted 
n e (cm"" 3 ) 

• 100 fj,m Emissivity 
per H + Ion a 

100 ‘/rm Emissivity 
per H atom a 

Reference 

HQBN 

48-145 

60-75 

0.08 

-6.1 ±6.8 

23.1 ±1.2 

This paper 

LH 

130-163 

46-65 

0.08 

8.9 ±5.9 

18.8 ±2.0 

This paper 

NEP 

93-100 

25-34 

0,08 

0.8 ±3.1 

18.9 ±2.0 

This paper 

Qi 

o 

CO 

i 

o 

00 

30-41 

0.08 

3.2 ±4.3 

20.8 ±1.7 

This paper 

Q2 (HI cut) 

97-170 

28-52 

0.08 

-1.2 ±2.9 

15.4 ±0.9 

This paper 

Q2 (no HI cut) 

97-170 

22-52 

0.08 

16.2 ±0.5 

14.9 ±0.1 

Lagache et al. (2000) 

l = 144, b = -21 

136-150 

16-27 

0.09 

-0.4 ±1.9 

16.3 ±1.1 

Arendt et al. (1998), this paper 

Spica 

297-343 

39-70 

0.6 

87 ±20 

19.5 ± 1.6 

Arendt et al. (1998), this paper 

Eridanus Superbubble 

180-210 

15-50 

0.8 

20 

19 

Heiles et al. (1999) 

Barnard’s Loop ■ 

208-217 

13-27 

2 

34 

17 

Heiles et al. (2000) 

Outer Galactic Plane 

90-270 

0-5 

2 

54 ±5 

6.8 ±0.4 

Sodroski et al. (1997) 

Inner Galactic Plane 

270-90 

0-5 

10 

380 ± 17 

37 - 180 

Sodroski et al. (1997) 


in nW m 2 sr 1 /10 20 cm -2 
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Fig. 1. — - Location of the regions analyzed on a Galactic coordinate Mollweide projection 
centered at l — 0°. The Lockman Hole (LH) region is the region mapped in H I by Snowden 
et al. (1994). The north ecliptic pole (NEP) region is the region mapped in H I b.y Elvis, 
Lockman and Fassnacht (1994). The DIRBE high quality B north (HQBN) region is defined 
as the region at Galactic latitude b > +60° and ecliptic latitude /3 > +45°. The first quadrant 
(Ql) region is defined by 30° < l < 80°, 30° < b < 41°. The second quadrant (Q2) region is 
a region previously studied by Lagache et al. (2000). 
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(d) lOOfj-m Intensity minus H I correlated component 



ao 70 60 60 40 so 

Galactic Longitude 


Fig. 2. — Images of the Ql region in (a) DIRBE 100 fan intensity after subtraction of 
the interplanetary dust emission model, (b) H I column density from the Leiden-Dwingeloo 
survey, (c) Ha intensity from the WHAM Northern Sky Survey, and (d) 100 fan intensity as 
in (a) but with the component that is correlated with H I column density subtracted. The 
distributions of H I and Ha differ from each other and differ from an isotropic distribution, 
so the infrared data can be decomposed into a sum of the three distributions. No correlation 
is seen between Ha and the residual 100 fim emission in (d), consistent with the low value 
of the Ha coefficient Hi obtained from our analysis. Possibly the 100 /mi emissivity of the 
ionized medium is low and the intensity variations in (d) are not related to the ionized 
medium, or Ha is not a good tracer of the 100 fan emission from this medium. The image 
display ranges, from black to white, are 20 to 170 nW m~ 2 sr -1 for (a), 0.4 to 6.3 10 20 atoms 
cm -2 for (b), 0.5 to 1.6 R for (c), and -25 to 50 nW m -2 sr -1 for (d). The black contour in 
(b) shows the 1V(H I) cut of 3 10 20 atoms cm -2 used in the analysis, as described in §3.1. 


/(Ha) (Rayleighs) ^(100 /tm) (nW 
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0 1 2 3 4 0.0 0.5 1.0 1.5 2.0 

AT(H I) (10 20 H atoms cm’ 2 ) /(Ha) (Rayleighs) 


Fig. 3. — Correlation plots for Ql region positions used in our analysis, below an N ( H I) cut 
at 3 10 20 atoms cm -2 . The correlation between 100 /nil intensity and H I column density 
(a) is tighter than that between 100 //in intensity and Ha intensity (b) or that between 
Ha intensity and H I column density (c). No correlation is seen between Ha intensity and 
residual 100 pm intensity after subtraction of the H I correlated component (d). The cross 
plotted in the lower right of each panel shows the typical statistical measurement uncertainty 
(± lcr) for each quantity. 
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Fig. 4. — Fits to the 100 // m data for the Ql and Q2 regions. For each region, the 100 fim 
intercept is nearly the same for the fit using H I (top) and the fit using H I and Ha- (bottom). 
The scatter about the fit line is also nearly the same for the two cases. 
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Fig. 5. — Infrared emissivity per H nucleus for the neutral atomic gas phase (filled symbols) 
and the ionized gas phase (open symbols) from the three-component fits for each region. 
For the ionized gas phase, the emissivity values were obtained from the B r parameter values 
using a conversion factor of I(Ha)/N(H + ) = 1.15 Rayleighs/ 10 20 cm -2 (see text). The 
emissivity values for the ionized phase are consistent with zero, and the emissivity values 
for the neutral phase are consistent with those derived when an Ha-correlated component is 
not included in the fits (see Table 3). 
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Fig. 6. — Mean residual infrared intensity C x as a function of wavelength, from the three- 
component fits for each region. The residual intensity values for the five regions are consistent 
with isotropy at 140 and 240 //m, marginally consistent with isotropy at 100 /un, and not 
consistent with isotropy at 60 /nn. The error bars show la statistical uncertainties and 
do not include systematic uncertainties that contribute to the total uncertainty for a CIB 
measurement. 
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Fig. 7.— Residual infrared intensity averaged over the HQBN, LH, NEP, Ql, and Q2 
regions as a function of wavelength. Results from the two-component fits (circles) and 
from the three-component fits (squares) are in close agreement. The error bars show la 
statistical uncertainties and do not include systematic uncertainties that contribute to the 
total uncertainty for a CIB measurement. 
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Fig. 8. — Derived values of 100 /im emissivity per H + ion for different regions plotted as a 
function of estimated electron density. The data are from Table 7. Most previously studied 
regions have emissivities greater than the 2a upper limits shown for our regions. In most 
cases this can be attributed either to enhanced dust heating by nearby stars or, for regions 
with electron density of about 1 cm -3 or greater, to enhanced dust heating by Ly a radiation. 



